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INTRODUCTION 


Laser  diagnostic  techniques  have  been  developed  at  Yale  University  that  are  capable 
of  two-dimensional  mapping  of  multiple  scalars  and  velocities  in  turbulent  flames.  The 
techniques  are  designed  to  measure  quantities  and  flow  configurations  of  current  interest 
to  combustion  modelers.  The  availability  of  quantitative  data  on  the  spatial  and  temporal 
characteristics  of  structures  in  turbulent  reacting  flows  will  aid  in  understanding  the 
interaction  of  chemical  reactions  with  the  turbulent  motion.  A  better  understanding  of  this 
key  interaction  is  important  for  testing  existing  models  of  turbulent  combustion  as  well  as 
for  suggesting  new  models. 

During  the  three  years  of  AFOSR  support,  significant  progress  has  been  made  in 
refining  the  mixture  fraction  imaging  technique  and  in  two  critical,  related  research  areas: 
(1)  Three-scalar  measurements  in  turbulent  flames,  and  (2)  Optical  flow  velocimetry 
techniques  based  on  unseeded  scalar  image  pairs.  In  the  following  sections,  the  major 
accomplishments  are  outlined. 


RESEARCH  ACCOMPLISHMENTS 


Three-Scalar/Mixture  Fraction  Imaging 

During  the  current  funding  period,  significant  advances  have  been  made  in  the 
refinement  of  mixture  fraction  imaging  techniques  in  turbulent  nonpremixed  flames.  The 
mixture  fraction  (£)  is  defined  as  the  mass  fraction  of  atoms  originating  in  the  fuel  stream 
and  is  independent  of  the  chemical  reaction  occurring  in  the  flame.  The  gradient  of  the 
mixture  fraction  is  needed  to  find  the  scalar  dissipation,  which  determines  the  rate  of 
molecular  mixing  in  the  flow.  The  mixture  fraction  and  its  gradient  are  important 
parameters  in  the  modeling  of  turbulent  flames,  and  the  experimental  determination  of 
mixture  fraction  over  a  wide  field  is  essential  for  testing  these  models. 

One  method  for  constructing  a  conserved  scalar  suitable  for  imaging  experiments  in 
reacting  flows  has  been  through  the  simultaneous  measurement  of  temperature  (T)  and 
fuel  concentration  [1,2].  The  conserved  scalar  /J,  is  defined  based  on  fuel  mass  fraction 
(Yf)  and  enthalpy  and  takes  the  form: 

P  =  YF+cpT/Q  (1) 


where  Q  is  the  lower  heat  of  combustion  and  cp  is  the  specific  heat  at  constant  pressure. 
This  can  be  cast  into  an  expression  for  mixture  fraction: 
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This  two-scalar  approach,  which  assumes  unity  Lewis  number  and  idealized  one-step 
reaction  between  fuel  and  oxidizer,  relates  to  the  measured  signals  through: 


£FT=C£_Rm  C, 

WRa  Q 


Ra 


Cp.air^air 


(3) 


where  Rm  is  the  measured  fuel  Raman  scattering  and  Ra  is  the  Rayleigh  scattering.  The 
parameter  <7 ,  which  is  proportional  to  the  Rayleigh  cross  section,  the  mixture  molecular 
weight  W,  and  the  specific  heat  cp,  are  dependent  on  the  mixture  fraction.  Strained 
counterflow  flame  calculations  provide  appropriate  functional  forms  for  these  parameters 
which  are  incorporated  into  an  iterative  scheme  for  determining  £.  The  remaining 
constants  C,  and  C2  must  be  determined  from  calibration  experiments. 

In  the  mixture  fraction  imaging  experiments,  the  fuel  concentration  has  been  obtained 
using  Raman  scattering  from  the  fuel.  The  use  of  Raman  scattering  to  determine  the  fuel 
concentration  allows  the  use  of  chemically  simple  fuels  such  as  methane  and  hydrogen, 
which  are  reasonably  consistent  with  the  one-step  chemistry  assumptions  implicit  in  the 
technique.  Earlier  experiments  using  fluorescence  from  more  complex  hydrocarbons 
provided  good  signal/noise,  but  suffered  inaccuracies  in  the  derived  mixture  fraction  due 
to  loss  of  parent  fuel  [3,4]. 

In  a  previous  funding  period,  experiments  were  performed  in  a  variety  of 
nonpremixed  flame  configurations,  including  a  lifted  methane  diffusion  flame  [5].  Data 
sets  were  also  obtained  for  piloted  methane/air  flames  [6]  and  an  unpiloted  hydrogen 
flame  [7].  In  the  methane/air  flames,  the  methane  was  diluted  with  air  (70%  by  volume) 
to  eliminate  soot,  and  resulted  in  a  stoichiometric  mixture  fraction  of  E,  =  0.29.  The 
premixed  annular  pilot  flame  was  a  stoichiometric  mixture  of  acetylene,  hydrogen,  and 
air,  which  insured  a  steadily  burning  main  flame  without  significant  local  extinction.  The 
fuel  jet  issued  from  a  nozzle  of  diameter  6.1  mm  into  a  filtered,  vertical,  coflowing  7.0 
m/s  air  stream. 


One  major  issue  resolved  within  this  funding  period  concerns  uncertainty  in  the 
position  of  the  stoichiometric  contour  using  the  two-scalar  approach  with  air-diluted 
methane  flames.  While  dilution  prevents  difficulties  arising  from  excessive  luminosity 
(or  laser-induced  incandescence),  it  also  creates  a  broader  minimum  for  the  Rayleigh 
scattered  signal  in  mixture  fraction  space  as  shown  in  Fig.  1.  Close  to  the  stoichiometric 
value  of  the  mixture  fraction,  the  fuel  concentration  approaches  zero  and  the  Rayleigh 
signal  remains  nearly  constant.  It  was  proposed  that  a  suitable  third  scalar  quantity  would 
clarify  the  accuracy  of  the  two-scalar  approach.  This  could  be  achieved  either  by 
modifying  the  two-scalar  approach  or  providing  simultaneous  independent  determination 
of  mixture  fraction.  A  detailed  study  of  candidate  third  scalars  led  to  the  in-depth 
examination  of  three  particular  scenarios:  (1)  OH  fluorescence,  (2)  seeded  sodium  atomic 
fluorescence,  and  (3)  nitrogen  Raman  with  a  customized  fuel  mixture. 
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Figure  1.  Profiles  of  signal  intensities:  Rayleigh  scattering  and 
methane  Raman  scattering  scattering  as  a  function  of 
mixture  fraction  for  an  air-diluted  (3:1  by  volume)  methane 
counterflow  flame.  The  signals  are  scaled  to  appear  on  the 
same  plot. 


OH  Fluorescence  Imaging 

The  first  case  employed  laser-induced  fluorescence  of  the  OH  radical  as  an 
approximate  marker  of  the  reaction  zone  near  stoichiometric.  Experiments  were 
performed  in  a  Re= 15,000  turbulent  nonpremixed  methane  flame.  The  dye  laser  output 
was  frequency  doubled  and  tuned  to  excite  the  Q,(7)  transition  of  OH  (-289  nm)  and 
fluorescence  was  observed  at  a  wavelength  near  309  nm.  Laser  energy  in  the  ultraviolet 
was  about  2  mJ  per  pulse.  Simultaneous  Rayleigh/Raman/OH  LIF  images  taken  25  jet 
diameters  downstream  are  shown  in  Fig.  2.  The  imaged  region  is  20.5  mm  x  4.4  mm, 
which  includes  the  jet  centerline  on  the  left  side  of  the  image  and  a  region  of  ambient  air 
on  the  other  side  of  the  image. 

In  general,  the  OH  fluorescence  does  not  fill  regions  where  the  Rayleigh  signal  is  at  a 
minimum  and  Raman  signal  disappears.  While  this  technique  offers  insight  into  the  local 
flame  structure,  it  fails  to  provide  a  clear  advantage  in  terms  of  supplementing  mixture 
fraction  information  that  can  be  calculated  from  the  Rayleigh  and  Raman  signals. 
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Figure  2.  (a)  Rayleigh,  (b)  Raman,  and  (c)  OH  fluorescence  images 

from  a  turbulent  nonpremixed  methane  flame  (Re= 15000) 
taken  25  jet  diameters  downstream. 


Sodium  Fluorescence  Imaging 

In  the  past,  sodium  has  been  added  to  flames  for  applications  such  as  temperature 
measurements  (by  the  line-reversal  technique)  and  control  of  soot  formation.  Laser- 
induced  fluorescence  of  sodium  is  well  documented  [8]  although  the  present  work 
represents  the  first  attempt  to  assess  its  applicability  as  a  marker  of  mixture  fraction 
through  the  flame  zone.  Sodium  is  introduced  into  the  flow  using  an  atomizer  (TSI  Model 
9306)  containing  a  3.4xl0"3  M  solution  of  NaCl  in  initially  de-ionized  water.  The  NaCl 
molecules  dissociate  due  to  high  temperatures  encountered  in  the  flame  zone.  The 
Nd:  YAG-pumped  dye  laser  is  tuned  to  the  32S1/2 -*32P1/2itransition  of  sodium  at  589  nm, 
with  output  beam  energies  around  15  mJ.  The  resulting  image  represents  fluorescence  of 
the  32P 1/23/2  sodium  doublet  to  the  ground  state  at  wavelengths  of  589  nm  and  589.6  nm, 
respectively. 

Measurements  were  again  made  in  a  turbulent  flame  (/?e=  15,000)  25  jet  diameters 
downstream  as  shown  in  Fig.  3.  To  a  greater  extent  than  with  OH  LIF,  the  sodium 
fluorescence  signal  fills  the  region  between  the  Raman  and  Rayleigh  signal,  and  extends 
into  the  fuel  region.  The  peak  of  the  fluorescence  occurs  in  the  high  temperature  region, 
probably  near  the  stoichiometric  contour  given  the  results  of  laminar  flame  validation 
studies.  However,  Mie  scattering  from  residual  salt  particles  interferes  significantly  with 
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the  fuel  Rayleigh  signal  prohibiting  quantitative  determination  of  mixture  fraction.  This 
interference  also  prevents  contour  smoothing  the  Raman  image. 


Centerline  Min  Max 

Figure  3.  a)  Rayleigh,  (b)  Raman,  and  (c)  sodium  fluorescence 
images  from  a  turbulent  nonpremixed  methane  flame 
(Re= 15000)  taken  25  jet  diameters  downstream. 


Nitrogen  Raman  Imaging 

The  ideal  third  scalar  is  a  monotonic  function  of  mixture  fraction,  a  trait  not  exhibited 
by  either  OH  or  sodium.  An  inert  species  would  have  this  property  if  it  was  present  only 
in  either  the  fuel  or  oxidizer  stream.  This  concept  led  to  a  comprehensive  study  of 
imaging  nitrogen  mass  fraction  as  a  third  scalar  measurement  using  Raman  scattering 
with  a  customized  fuel  mixture  [9].  Experiments  were  performed  in  which  the  fuel  stream 
contained  no  nitrogen,  which  gave  the  nitrogen  signal  the  desired  characteristic  as  shown 
in  Fig.  4. 

Assuming  no  significant  nitrogen  consumption  occurs  during  reaction,  we  can  write 
the  conserved  scalar  in  terms  of  nitrogen  mass  fraction: 

P"'  -  YKi  (4) 
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with  the  mixture  fraction: 


^  =1- 


'2  _ 


v 

1N2,aii 


=  1- 


c3  CT 


(5) 


where  RmNi  is  the  nitrogen  Raman  signal,  and  C3  is  a  calibration  constant.  This 

formulation  requires  measurement  of  the  temperature  and  nitrogen  concentration,  and 
thus,  represents  an  additional  two-scalar  approach. 


Mixture  Fraction 

Figure  4.  Profiles  of  signal  intensities:  Rayleigh  scattering,  methane 
Raman  scattering,  and  nitrogen  Raman  scattering  as  a 
function  of  mixture  fraction  for  an  argon/oxygen-diluted 
methane  counterflow  flame.  The  signals  are  scaled  to 
appear  on  the  same  plot. 


The  customized  fuel  used  in  the  experiments  consists  of  methane  mixed  with  argon 
and  oxygen  to  provide  an  overall  3/1  dilution  ratio  by  volume  (3  parts  diluent  to  1  part 
fuel).  The  oxygen  content  was  selected  to  simulate  air  dilution,  which  gives  a 
composition  of  25%  methane,  59%  argon,  and  16%  oxygen  by  volume.  This  fuel  mixture 
has  a  stoichiometric  mixture  fraction  4=0.41,  putting  the  reaction  zone  well  inside  the 
shear  layer.  Counterflow  flame  calculations  (strain  rate=100  s'1)  of  this  system 


6 


demonstrate  the  departure  of  the  ^  two-scalar  formulation  from  the  “actual”  mixture 
fraction  (Fig.  5)  calculated  using  the  following  formula  [10]: 

2ZcIVc+iZKIVH+(Z0jlr-Z0)IW<> 

2Zc,m/Wc  +  iZ +(ZMr  -  Z0ilJl/W0 

It  is  predicted  from  these  calculations  that  the  mixture  fraction  based  upon  nitrogen 
mass  fraction  should  exhibit  little  deviation  from  the  actual  Deviation  from  one-step 
chemistry  (i.e.  loss  of  parent  fuel  to  intermediate  species)  is  compensated  in  calculating 
fuel  concentration  by  the  expression  RmF  =  RmF0  (1  +  CF  <I>2),  where  RmFfi  is  the 
uncompensated  value,  and  the  reactivity  is  defined  as  <E>  =  (T  -  TambientV(Tad 

”  ^ambient)  [61 

CF  is  a  weighting  coefficient  (CF  =  0.8)  and,  as  illustrated  in  the  figure  by  the  curve 
marked  corrects  on  the  rich  side  of  stoichiometric. 


Actual  Mixture  Fraction 

Figure  5.  Mixture  fraction  calculated  from  strained  laminar  flame 

calculations  (100s1)  using  the  fuel-temperature  (black  line) 
and  nitrogen-temperature  (blue  line)  two-scalar  approaches 
plotted  against  mixture  fraction  calculated  using  Equation 
3.  The  effect  of  fuel  correction  on  £FT  is  also  shown  (red 
line). 


Simultaneous  planar  Rayleigh,  fuel  Raman,  and  nitrogen  Raman  images  have  been 
collected  in  experiments  using  three  cameras  and  a  single  laser.  Figure  6  shows  the 
experimental  facility.  A  flashlamp-pumped  dye  laser  is  employed  in  an  intracavity 
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configuration  to  generate  single-shot  energies  up  to  4.7  J  at  532  nm.  The  beam  is  focused 
into  a  sheet  over  a  6.1  mm  diameter  piloted  burner.  Scattered  light  is  collected  on  both 
sides  of  the  flame  by  low  //#  camera  lenses  oriented  perpendicular  to  the  laser  sheet.  The 
Rayleigh  scattering  and  fuel  Raman  scattering  are  collected  along  the  same  optical  path 
and  divided  with  a  50/50  pellicle  beam  splitter,  while  the  weaker  nitrogen  Raman 
scattering  is  collected  along  the  opposite  optical  path.  Image  intensifiers  are  lens-coupled 
to  liquid-cooled  CCD  cameras  and  isolated  with  appropriate  interference  filters. 


Sheet-Forming  Lens 


Flashlamp  Pumped  Dye  Laser 
(532  nm) 


N2  Raman 


CCD  Camera 


// 1.4  Lens 
Beam  Splitter 


Burner 


image  Intensifier 
Interference  Filter 
//1 .2  Lens 

Cylindrical  End  Mirror 


Rayleigh 


CH4  Raman 


CCD  Cameras 


Figure  6.  Schematic  of  the  three-scalar  turbulent  flame  imaging 
experiment. 


In  the  past,  the  imaging  resolution  of  the  optical  setup  has  been  characterized  in  terms 
of  the  volume  associated  with  each  pixel.  For  the  present  work,  the  pixel  volume  is  49  x 
49  x  500  pm3,  where  the  largest  value  corresponds  to  the  laser  sheet  thickness,  and  the 
remaining  dimensions  describe  the  area  imaged  onto  a  single  pixel.  The  actual  spatial 
resolution  is  a  more  complex  function  of  the  optical  layout,  including  alignment,  lenses, 
filters,  image  intensifiers,  and  camera  pixel  size.  In  order  to  better  quantify  the  spatial 
resolution  for  this  configuration,  simultaneous  images  have  been  taken  of  a  uniformly 
illuminated  25  pm  wire  located  at  the  focal  plane.  Individual  camera  resolutions  are 
based  upon  the  resulting  full-width  half-maximum  (FWHM)  intensity  of  the  wire  image 
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from  each  camera  following  scaling,  translation,  rotation,  and  cropping.  The  spatial 
resolution  is  170  pm  on  the  Rayleigh  camera,  140  pm  on  the  nitrogen  Raman  camera, 
and  275  pm  on  the  fuel  Raman  camera.  Translating  the  wire  normal  to  the  plane  of  the 
laser  sheet  within  the  beam  thickness  (±  250  pm  from  the  focal  plane)  has  minimal  effect 
(<10%)  on  the  individual  camera  resolutions. 

Because  this  is  a  multi-camera  experiment  and  the  cameras  are  located  along  different 
optical  trains,  there  remains  the  issue  of  how  well  the  images  correlate  on  a  pixel-by- 
pixel  basis.  Cross-camera  spatial  resolutions  here  are  defined  based  upon  adding 
combinations  of  matched  images  from  each  camera,  and  measuring  the  resultant  full- 
width  half-maximum  of  the  wire.  With  optimal  matching,  the  on-axis  spatial  resolution  is 
280  pm,  about  equal  to  the  largest  single  camera  resolution.  Factors  such  as  distortion 
may  cause  degradation  away  from  the  optical  axis  to  a  maximum  measured  400  pm. 
Using  the  He2  intensity  point  rather  than  the  FWHM  for  determining  spatial  resolution 
increases  the  reported  values  by  -75%. 

Estimates  of  the  Kolmogorov  scale  (k)  on  the  centerline  for  a  diluted  methane  flame 
(3/1  air/methane  by  volume)  give  a  value  of  k=95  pm  at  /te=20,600  [4].  Using  this  as  an 
approximate  value  for  the  turbulent  flame  in  this  work,  the  resolutions  reported  here  are 
in  the  range  2-4k.  A  study  of  an  isothermal  jet  performed  by  Namazian,  et  al.  reports  that 
a  spatial  resolution  of  5k  should  be  sufficient  for  capturing  60%  of  the  scalar  dissipation 
spectrum  [11],  In  flames,  where  heat  release  is  expected  to  increase  length  scales,  a 
resolution  of  5k  should  be  sufficient  to  record  most  of  the  scalar  dissipation. 

Figure  7  shows  line  plots  of  mixture  fraction  from  a  laminar  flame  (Re=1600)  at  a 
location  15  nozzle  diameters  downstream  (D=6.1  mm).  Two  curves  are  shown  for  £FT, 
which  differ  in  the  parameterization  of  the  mixture  fraction  dependent  terms  appearing  in 
Eqn.  3  (i.e.  o,  W,  cp).  The  curve  marked  “No  lean  correction”  uses  flame  calculation 
terms  parameterized  by  “actual”  mixture  fraction  determined  from  the  Bilger  formula 
(Equation  6). 

To  correct  £Fr  on  the  lean  side,  cp(£),  a(£),  and  W(£)  are  described  as  functions  of  the 
predicted  two-scalar  mixture  fraction  ^FT,  rather  than  the  actual  mixture  fraction  from  the 
flame  calculations.  Figure  5  showed  that  the  old  approach  results  in  a  departure  of 
from  the  more  rigorous  formulation,  as  indicated  by  strained  laminar  flame  calculations 
(100  s'1  strain  rate).  Deviation  from  one-step  chemistry  (i.e.  loss  of  parent  fuel  to 
intermediate  species)  is  compensated  in  the  curve  marked  by  using  a  weighting  term 
involving  reactivity.  The  curve  for  £Nz  exhibits  little  deviation  from  the  actual  mixture 
fraction.  By  assigning  functional  dependences  based  on  the  predicted  £FT  from  flame 
calculations,  the  t**  curve  shown  in  Fig.  7  is  obtained  showing  improved  agreement  with 

The  flame  computations  are  insensitive  to  variations  in  the  strain  rate  over  the  range 
10-200  s'1,  which  is  expected  to  be  representative  of  the  scalar  dissipation  values 
measured  in  the  turbulent  flame,  based  on  existing  data  in  similar  flames  [6,12]. 
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Figure  7.  Radial  variation  of  mixture  fraction  from  the  modified  fuel- 
temperature  (black  line)  and  nitrogen-temperature  (blue 
line)  two-scalar  approaches  in  a  laminar  flame.  The 
unmodified  approach  overpredicts  mixture  fraction  in 
regions  around  stoichiometric  (red  line). 


Applying  this  technique  to  single-pulse  imaging  in  turbulent  nonpremixed  flames 
provides  similar  results,  although  the  N2  Raman  signal  is  more  affected  by  noise, 
especially  in  regions  of  high  mixture  fraction  where  there  is  little  nitrogen.  Figures  8  and 
9  show  images  taken  25 D  downstream  from  a  Re= 15,000  turbulent  flame;  the  Raman 
images  have  been  contour  smoothed  [13].  Qualitatively,  the  scalar  dissipation  fields,  % 

(defined  as  x  =  2©  V£  •  V£,  where  27is  the  diffusivity),  appear  similar,  revealing  the 
same  main  structural  features,  with  significant  scalar  dissipation  apparent  along  the  edge 
of  the  main  jet.  The  position  of  the  stoichiometric  mixture  fraction  contour  is  highlighted 
(white  lines)  in  Fig.  9. 

The  experiments  measuring  nitrogen-temperature  mixture  fraction  have  served  as  a 
guide  for  correcting  the  fuel-temperature  mixture  fraction  for  values  around  and  lean  of 
stoichiometric.  It  has  been  shown  that  parameterizing  specific  heat,  molecular  weight, 
and  Rayleigh  cross  section  as  a  function  of  ^  predicted  from  counterflow  flame 
calculations,  rather  than  the  actual  mixture  fraction,  improves  the  performance  of  this 
two-scalar  approach.  Under  turbulent  conditions  (Re=15,000),  the  two  approaches  reveal 
differences  close  to  the  centerline,  most  likely  a  result  of  noise  limitations  of  the  nitrogen 
Raman  signal.  This  work  increases  confidence  in  employing  ^  for  mixture  fraction 
determination,  which  remains  the  most  attractive  approach  because  of  its  superior  signal- 
to-noise  characteristics. 
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Instantaneous  Rayleigh,  CH4  Raman,  and  N2  Raman  images 
of  the  turbulent  flame. 
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Optical  Flow  Velocimetry  ( OFV) 

Extending  the  mixture  fraction  imaging  approach  to  include  simultaneous  planar 
velocity  data  is  extremely  desirable.  Such  measurements  would  give  direct  information 
on  scalar  fluxes  and  would  be  valuable  for  validating  model  assumptions.  The  availability 
of  multi-dimensional  velocity  and  scalar  measurements  would  allow  further 
characterization  of  the  flames  by  providing  information  on  the  relationship  of  scalars  and 
their  gradients  with  the  vorticity  and  strain  rate. 

Planar  imaging  of  the  velocity  field  in  reacting  and  non-reacting  flows  is  typically 
obtained  by  imaging  of  Mie-scattering  from  particles,  either  with  PIV  or  particle  tracking 
velocimetry  (PTV).  In  a  previous  funding  period,  we  demonstrated  that  PIV  can  be 
combined  with  scalar  measurements  using  PLIF  [14].  However,  particle-based  techniques 
complicate  simultaneous  laser-based  measurements  of  certain  scalar  quantities  due  to 
interference  from  the  seed  particles.  In  particular,  simultaneous  imaging  of  Rayleigh 
scattering  is  impossible,  making  temperature  and  mixture  fraction  measurements 
difficult.  A  number  of  “non-particle”  velocimetry  approaches  have  been  applied  in 
turbulent  flows,  though  few  have  found  broad  applicability  under  reacting  conditions  [15- 
17]. 

Optical  flow  is  concerned  with  the  apparent  motion  of  an  image  (or  brightness) 
pattern  between  two  frames  [18].  Under  well-posed  conditions,  the  optical  flow  will 
faithfully  characterize  the  actual  motion  field  giving  rise  to  the  sequential  images.  Several 
optical  flow  methodologies  exist,  and  the  more  common  approaches  are  based  on  either 
differential  techniques  or  region-based  matching  [19].  In  recent  work,  we  have  examined 
region-based  approaches  [20]  because  these  do  not  suffer  from  the  difficulties  associated 
with  numerical  differentiation  or  the  requirement  of  small  displacements  (~1  pixel/frame) 
as  is  the  case  with  differential  techniques.  These  and  other  approaches  can  at  best  provide 
an  approximation  to  the  instantaneous  two-dimensional  fluid  velocity  given  a  planar 
scalar  measurement.  The  approximation  is  valid  when  the  velocity  in  the  third, 
unmeasured  direction  is  small,  or  the  scalar  gradient  in  that  dimension  is  small  [21]. 

Although  the  various  approaches  are  sometimes  calibrated  using  synthetic  images 
[22],  direct  comparisons  between  velocities  obtained  from  optical  flow  and  PIV  are  not 
available.  In  an  attempt  to  fill  this  gap,  we  performed  preliminary  experiments  in  an 
isothermal  turbulent  jet  of  air  seeded  with  both  acetone  and  alumina  particles.  The 
experimental  configuration  is  shown  in  Fig.  10.  The  nozzle  used  in  these  experiments 
consists  of  simple  tube  with  a  squared-off  end  (6.35  mm  OD;  4.57  mm  ID).  The  turbulent 
jet  issues  into  an  unconfined  low  velocity  coflow  of  air  (~1  m/s).  The  illumination  source 
consists  of  the  fourth  harmonic  (266  nm)  of  an  Nd:YAG  laser  (Continuum,  Powerlite 
8000)  which  is  formed  into  a  sheet  approximately  5  mm  high  with  a  beam  energy  of  ~2.5 
mJ  per  pulse.  The  laser  is  operated  in  dual-pulse  mode  with  a  pulse  separation  of  57  |xs. 
The  camera  employed  for  the  acetone  fluorescence  is  a  dual-frame  Cooke  Corp./PCO 
SensiCam  with  an  interline  progressive-scan  CCD  (1280  x  1024  pixels).  This  camera 
allows  the  acquisition  of  two  exposures  separated  by  as  little  as  200  ns.  The  primary 
collection  lens  is  a  large-format  camera  objective  (Canon,  85  mm,  //1.2)  with  a 
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corresponding  magnification  of  31  pixels/mm.  A  glass  filter  (Schott,  WG305)  is  used  to 
minimize  residual  Mie  scattering.  The  acetone  fluorescence  is  broadband,  peaking  around 
430  nm.  This  setup  allows  the  acquisition  of  two  acetone  fluorescence  images  for  use  in 
optical  flow  velocimetry  (OFV).  A  Photometries  CH350  camera  with  a  UV  sensitive 
CCD  coating  (QE=10%)  and  a  Nikon  Nikkor  UV  lens  (//4.5)  allows  direct  imaging  of 
the  Mie  scattering  from  the  seed  particles.  The  two  particle  scattering  frames  are  spatially 
separated  on  the  CCD  by  operating  a  rotating  mirror  at  approximately  50  Hz. 


Cross¬ 

correlation 

PIV 


DC  Motor 


Figure  10.  Experimental  configuration  for  the  two-camera  cross¬ 
correlation  PIV  and  acetone  fluorescence  scalar  imaging 
OF  velocimetry  validation. 


Single-shot  scalar  images  from  the  turbulent  non-reacting  flow  (/?e=2810)  are  shown 
in  Fig.  11,  where  the  color  scale  indicates  concentration.  From  this  pair  of  images,  the 
corresponding  PIV  results  and  OFV  using  two  different  algorithms  are  shown  in  Fig.  12. 
Further  details  regarding  the  algorithms  will  be  available  in  an  upcoming  paper  [23]. 
Qualitatively,  the  results  are  promising,  however,  further  detailed  analysis  is  required  to 
determine  the  strengths  of  particular  algorithms  and  lead  to  a  unified  OFV  approach 
suitable  for  a  wide  range  of  conditions  including  reacting  flows. 
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Figure  11.  Scalar  fluorescence  images  taken  in  an  turbulent, 
isothermal  jet  of  acetone-seeded  air. 


In  summary,  proof-of-concept  experiments  were  performed  within  the  current 
funding  period  to  compare  simultaneous  optical  flow  velocimetry  (OFV)  and 
conventional  PIV.  In  OFV,  velocity  vectors  are  computed  based  upon  the  displacement 
fields  of  molecular  flow  tracers  which  would  allow  concurrent  measurement  of  mixture 
fraction  using  the  two-scalar  approach.  Significant  progress  has  been  made  in 
implementing  exemplar  cross-correlation  based  OFV  algorithms,  and  future  work  in  the 
areas  of  diagnostics  and  computational  methodology  will  result  in  rapid  advances  in  the 
field.  Current  results  have  been  made  available  to  other  researchers  via  the  World  Wide 
Web  to  have  broader  impact  and  encourage  new  development  [24]. 


t  6.0  m/s 


i  2  mm  j 


Figure  12.  Velocity  fields  corresponding  to  the  images  of  Fig.  1 1 .  Two 
different  optical  flow  velocity  fields  are  shown,  labeled 
Komiyama-OF  and  Fielding-OF  to  designate  the  algorithm 
employed.  The  third  vector  field  is  determined  using  PIV 
from  simultaneous  particle  images  (not  shown). 
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Additional  Research  Developments 

The  dual-frame  camera  used  in  the  OFV/PIV  experiment  offers  additional  capability 
beyond  velocimetry  applications  (both  PIV  and  optical  flow  techniques).  The  fast 
electronic  shuttering  (down  to  100  ns)  offers  the  possibility  of  eliminating  additional 
optical  components  necessary  in  intensified  systems  where  gateability  is  the  only 
requirement.  This  improves  the  spatial  resolution  of  the  imaging  system. 

To  demonstrate  this,  a  turbulent  air-diluted  methane  flame  (54%  CH4, 46%  air  by 
volume)  was  examined  with  the  non-intensified  imaging  system.  This  blue,  moderately 
sooty  flame  was  selected  in  order  to  assess  the  ability  to  discriminate  against  flame 
interference  using  the  camera’s  fast  shutter  capability.  For  this  experiment,  a  flashlamp- 
pumped  dye  laser  (Candela,  LFDL-20, 1.2  J/pulse)  was  employed  for  Rayleigh  scattering 
measurements.  The  luminosity  was  integrated  single-shot  for  4  ps  (to  bracket  the  duration 
of  the  2  ps  dye  laser  pulse)  and  1000  ps,  where  the  latter  value  was  intended  to  be 
representative  of  a  CCD  camera  operating  with  a  conventional  mechanical  shutter  and 
minimum  1  ms  exposure  time.  In  the  fast  shutter  case,  luminosity  accounts  for  less  than 
5%  of  the  signal  relative  to  a  standard  helium  background  image,  and  about  3%  of  the 
signal  in  the  hot  zone  of  the  flame.  This  contrasts  with  the  long  exposure  case  where 
these  values  increase  to  48%  and  25%,  respectively. 

Subsequent  studies  of  the  electronic  gating  capability  were  carried  out  during  laser- 
induced  incandescence  (LII)  experiments  in  a  pure  ethylene  laminar  flame.  The  LII  signal 
has  a  duration  of  several  hundred  nanoseconds,  and  an  intensifier  gate  time  of  ~1  ps  is 
common.  In  this  case,  the  extremely  bright  flame  luminosity  leaks  around  the  mask  of  the 
interline  CCD  frame,  increasing  the  background  contribution  to  the  imaged  signal  and 
consequently,  degrading  the  dynamic  range  of  the  weaker  LII  signal.  Thus,  current 
electronic  gating  features  of  fast-frame  cameras  must  be  employed  with  caution 
according  to  the  particular  application.  Unintensifed  Rayleigh  measurements  in  turbulent 
nonpremixed  flames  are  possible  when  using  the  common  air-dilution  to  minimize  soot 
production.  The  improvement  in  spatial  resolution  will  positively  impact  the  scalar 
dissipation  measurements  described  in  this  report. 
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